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ABSTRACT 

The Swift era has posed a challenge to the standard blast-wave model of 
Gamma Ray Burst (GRB) afterglows. The key observational features expected 
within the model are rarely observed, such as the achromatic steepening ('jet- 
break') of the light curves. The observed afterglow light curves showcase ad- 
ditional complex features requiring modifications within the standard model . 
Here we present optical/iVXff observations, millimeter upper limits and com- 
prehensive broadband modelling of the afterglow of the bright GRB 0505025A, 
detected by Swift. This afterglow cannot be explained by the simplistic form of 
the standard blast-wave model. We attempt modelling the multi-wavelength 
light curves using (i) a forward-reverse shock model, (ii) a two-component 
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outflow model and (iii) blast-wave model with a wind termination shock. The 
forward-reverse shock model cannot explain the evolution of the afterglow. 
The two component model is able to explain the average behaviour of the 
afterglow very well but cannot reproduce the fluctuations in the early X-ray 
light curve. The wind termination shock model reproduces the early light 
curves well but deviates from the global behaviour of the late-time afterglow. 

Key words: gamma-ray burst: individual : GRB 050525A 



1 INTRODUCTION 

Gamma Ray Bursts (GRBs) are extremely energetic cosmic explosions which outshine the 
entire 7-ray sky for a few seconds. The launch of Swift (Gehrels et al. 2004), a dedicated 
satellite to detect GRBs and rapidly follow-up their afterglow emission, has revolutionized 
the study of the most energetic cosmic explosions in the Universe. 

In the standard blast-wave model for GRB afterglows ( Rees fc Meszaros||1992 Paczynski 



& Rhoads 1993 also see Piran 1999 for a review), a relativistic shock decelerates through 



uniform circumburst medium, heats up the matter, accelerates particles and enhances the 
magnetic field downstream. Synchrotron radiation from the shocked particles is observed as 
the afterglow. The snap-shot synchrotron spectrum can be characterized by four spectral 
parameters (apart from the electron index p): the injection frequency z/ m , the cooling fre- 
quency z/ c , the self-absorption frequency v a , and the peak flux f m . The spectral parameters 
can be mapped to four physical parameters: the isotropic equivalent energy Ei SO , the ambient 
medium density (parameterized as number density no for a constant density medium and 
as A+ for a wind driven medium for which p(r) = 5.5 x 10 11 A * , (—) 

x ~\ / gem -1 Vcm/ 



as m 



Chevalier & 



Li 1999), and the fractional energy content in non-thermal electrons and magnetic field [e e 



and €b respectively). Jet-break, a simultaneous steepening seen in the multi-frequency light 



curves considered as a signature of the collimated outflow from the burst (Rhoads 1999), if 



observed (at time tj since burst), gives a handle on the initial collimation angle (9j) of the 
explosion, and hence to the total kinetic energy involved (-Etot)- 

The model was largely successful in explaining the pre-Swift observations of GRB after- 
glows. However, Swift with its ability to locate the afterglow within minutes of the burst and 
follow it up in UV, optical and X-ray bands has revealed complexity in the early afterglow 
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Multi-wavelength modelling of GRB 050525A afterglow 3 
emission that is not predicted by the model. The X-ray light curves in the Swift era have 



been rather dramatic, with steep decays, plateaus and flares, not witnessed earlier (Nousek 



et al. 2006 Chincarini et al. 2007). Moreover, in several afterglows, the flux evolution did 



not follow the predicted spectral-temporal relations (Liang et al. 2008). This has led to the 
conclusion that afterglow light curves differ drastically from burst to burst , owing to various 



physical processes shaping the flux evolution in various bands (Zhang et al. 2006). Another 
open issue in the Swift era is the absence of a jet-break. These complications often make it 
a demanding task to extract the physics of the burst and its surroundings from afterglow 
data. 



The bright low redshift (z = 0.606, Foley et al. 2005) Gamma Ray Burst 050525A was 



detected by the Swift-B AT on 2005 May 25 at 00:02:53 UT ( [Band et al.||2005| ). We refer 
to the burst trigger time as to- An isotropic equivalent 7-ray energy of 2.3 x 10 52 erg is 
inferred for the observed BAT fluence at a distance of 3.57 Gpc (assuming Q m = 0.3, Qa = 
0.7 and Hq = 70 km s _1 Mpc -1 ). The UVOT F-band observations started at T = t — 
t ~65 s and XRT observations began T ~ 75 s leading to well sampled early afterglow 
light curves. The proximity of the burst and the brightness of the afterglow made it a very 
good target for multi-wavelength observations. Ground based optical observations including 



robotic telescopes (Klotz et al. 2005 Delia Valle et al. 2006), radio observations in multiple 



frequencies by the Very Large Array (Cameron & Frail 2005) and Spitzer observations at 



~ 2 days in multiple _7?-bands (Heng et al. 2008) have been reported in the literature. 
However, a detailed modelling involving the full evolution of the relativistic shock to infer 
the physical parameters has not been attempted for this burst. 

In this paper we present a new set of VRIJH observations using eight different optical 
telescopes and millimeter upper limits from IRAM Plateau de Bure Interferometer (PdBI). 
We supplement our data with Swift-XJVQT and XRT data reported by Blustin et al.| ( p006 ), 
the optical data reported by Klotz et ab| ( |2005 ) and Delia Valle et al. (2006) and the radio 
data from the VLA afterglow repositorjj^to study the broad band evolution of the afterglow. 
We model the afterglow using three different extensions of the standard blast-wave model 



e.g. the forward-reverse shock model (Meszaros & Rees 1999), the two-component model 
( Berger et al.||2003 ) and a model including a wind termination shock ( Pe'er fc Wijers| [2006). 
Section [2] gives a description of the data acquired from different telescopes and the analysis 



1 http:/ /www. aoc.nrao.edu/~dfrail/allgrb_tablc.shtml 
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techniques. Multi-wavelength modelling under various premises are described in Section [3} 
Section [4] provides a summary of the multi-wavelength modelling results. 



2 DATA ACQUISITION AND REDUCTION 
2.1 Millimeter wave observations 

The IRAM Plateau de Bure Interferometer ( Guilloteau et al.]|1992 ) (PdBI, France) observed 
the position of GRB 050525A simultaneously at 92.682 GHz and 214.712 GHz about one 
day after the burst (see Table [J . The fields of view of the 15-m antennas at these frequen- 
cies are respectively 54'.'3 and 23'.'5, and the synthesized beams in the compact 5-antenna 
configuration were 7'.'05 x 5"58 at PA = 74° and 3"10 x 2/36 at PA = 78° where the position 
angle PA of the major beam axis is defined from North over East. The data reduction was 
done with the GILDAS^] software. We performed point source fits in the UV plane; for a 
fit fixed to the phase center position, a signal of at least 3a is needed to claim a detection, 
whereas a point source elsewhere in the field of view must have at least 5cr. We did not 
detect a millimeter counterpart, neither on the phase center coordinates nor elsewhere in 
the field of view. In the following, we take the 3cr levels as upper limits, i.e. 1.02 mJy for 
92.682 GHz and 3.66 mJy for 214.712 GHz (Table Q. 



2.2 Optical and NIR observations of the afterglow 

The afterglow of GRB 050525A was observed using different optical facilities in the broad 
bands VRLJH during 2005 May 25 to July 01. We used the 0.2m BOOTES1-B telescope 
located at the INTA-CEDEA station at El Arenosillo, Huelva, Spain; the 1.04m Sampur- 
nanand Telescope (ST) in India; the 1.2m semi-robotic Mercator telescope at La Palma; the 
1.5m telescope at the Observatorio de Sierra Nevada (OSN); the 2.01m Himalayan Chandra 
Telescope (HCT) in India; the 2.5m Issac Newton Telescope (INT) and the 2.2m and 3.5m 
telescopes at the Calar Alto Observatory (CAHA) for the observations (Table [2]). Several 
bias and twilight flat frames were acquired from different telescopes for pre-processing the 
CCD images. The pre-processing was done in a standard fashion including bias subtraction, 
flat fielding and cosmic ray removal in all object frames. Standard data reduction software 



2 http: / /www. iram.fr/IRAMFR/GILDAS 
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(IRA^\ MIDAlQ and DA OP HO 1f\ ( |Stetson|| 198 7[ ) ) were used for photometric analysis. The 
instrumental magnitudes of the optical afterglow were differentially calibrated using nearby 



secondary stars from the list of Henden (2005). A full compilation of VRIJH magnitudes of 
the afterglow is presented in Table [2} 

The earliest observations were taken with the BOOTES1-B telescope at T ~ 0.0044 days. 
We detect the afterglow in a co-added early time BOOTES1-B image. Figure [T] shows the 
optical afterglow in the early time BOOTES1-B image and in an image taken ~ 2 days after 
the burst with the OSN telescope. Late-time observations of the afterglow on 2005 June 
30 and 2005 July 01 with the 2.2-m CAHA and the 2.5-m INT respectively resulted in a 
non-detection and an upper limit. 

In Figure [2] we show the optical/NIR and mm data presented in this paper along with 
the published data in X-ray, optical and radio frequencies. We present an /-band light curve 
for this burst for the first time, which extends from T ~ 0.05 to 0.6 days and our data 
fills the gap in F-band light curve between 0.05 - 0.2 days (T ~ 5800 - 17,000 s). Our 
infrared observations (J- & i/-bands) are clustered around 0.09 days, but provide spectral 
information within the optical band that constrains model parameters. 



2.3 Host galaxy observations 

We obtained deep imaging of the GRB 050525A field in the gri bands between 2011 May 
26-28 with the 2.2m CAHA telescope. The data reduction was carried out in a manner 



similar to the one described in Section |2.2| A very faint source is detected in the i band at 
the location of the afterglow but only upper limits can be imposed in the g and r bands. 
The log of the observations is given in Table [2j We attempted further deep imaging in the J 
band on 2011 July 18-19 for a total exposure time of 2.33 hours spread over two consecutive 
nights with the 3.5m CAHA telescope equipped with the near-IR camera Omega2ooo- The 
photometric calibration is carried out by observing the UKIRT faint standards FS28 and 
FS30. We detect four faint sources close to the Swift XRT error circle in the J-band image. 

3 IRAF stands for Image Reduction and Analysis Facility distributed by the National Optical Astronomy Observatories which 
is operated by the Association of Universities for research in Astronomy, Inc. under co-operative agreement with the National 
Science Foundation 

4 MIDAS stands for Munich Image and Data Analysis System designed and developed by the European Southern Observatory 
ESO in Munich, Germany 

5 DA OPHOT stands for Dominion Astrophysical Observatory Photometry 
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The source closest to the enhanced XRT positions and UVOT position reported in Blustin 



et al. (2006) is indicated in Figure pi and is assumed to be the host galaxy of GRB 050525A. 



The host galaxy photometry is presented in Table [2j 



3 MULTI- WAVELENGTH MODELLING RESULTS 

GRB 050525A has a relatively well sampled afterglow light curve (see Figure [2]). The Swift 
observations end at T ~ 1 day while ground based i?-band observations continue up to 
~ 20 days. The early X-ray light curve starts off with a powerlaw decay but a minor flare 



is observed at T ~ 0.0035 days (Blustin et al. 2006). This flare is not seen in the optical 
light curves that undergo a smooth powerlaw decay. However, around T ~ 0.02 days, a 
re-brightening is observed in the i?-band light curve. There is considerable scatter in the 
i?-band light curve at this epoch and unfortunately no other frequencies cover this epoch. 
The afterglow evolution after this period appears to be within the expectations of the stan- 
dard blast-wave model with a near simultaneous break and similar post break indices in 



the optical and X-ray light curves (Blustin et al. 2006). However, the light curve before 



0.02 days significantly deviates from the extrapolation of the powerlaw decay that is present 
afterwards. There is a change in the i?-band light curve normalization after this epoch, or in 
other words, an increase in optical to X-ray spectral index with no color evolution within the 
X-ray or optical bands. This suggests that the value of v c decreases rapidly in this period. 
A single component, like a forward shock emission, alone cannot explain the complete evo- 
lution of this afterglow. Either there are multiple emission components or there is a density 
discontinuity in the ambient medium. 

We perform simultaneous multi-wavelength fitting of the afterglow data to obtain the 



underlying physical parameters. Parameterized extinction laws for Milky Way (Cardelli et al. 



1989), LMC, SMC (Pei 1992) and starburst galaxy (Calzetti 1997) are used for modelling 



extinction due to dust column in the host galaxy. A chi-square minimization procedure 
is used to derive the best fit parameters. For fitting we use all data except the late-time 
(T > 5 days) data in i?-band that is possibly contaminated by a supernova (SN) component 
(see section 



3.4) 



For the forward-reverse shock model and the two-component model, our code calculates 
the synchrotron flux f v (t) at an observed frequency v at a given time t as a function of 



http: / /www. swift, ac.uk/xrt_positions/ 
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the four spectral parameters which evolve following the respective shock dynamics (Wijers 



& Galama 1999; Resmi & Bhattacharya 2008). We assume both wind driven and constant 



density ambient medium while searching for the best fit. The free parameters of forward- 
reverse shock model and two-component model are the input values of u m , u c , v a and f m at 
a given epoch, the jet-break time tj, electron index p and the dust extinction E(B — V) in 
the host galaxy. If the electron distribution is a hard powerlaw, the synchrotron frequency 
Vi corresponding to the upper cut-off of the distribution and an index q parameterizing the 



time evolution of the cut-off will also be included as fit parameters (Resmi & Bhattacharya 



2008). In the forward-reverse shock model, the spectral parameters of the reverse shock 



and the shock crossing time t x are also included as free parameters. In the two-component 
model we have the deceleration time tdec of the wide jet as an additional free parameter. A 
band type smoothening is used in spectral (between multiple powerlaws of the synchrotron 
spectrum) and in the temporal (for gradual transition of the spectral parameters across tj, 



tdec and t x ) domains (Granot & Sari 2002). We derive the physical parameters from the 
best fit spectral parameters. For the wind termination shock model, we use the afterglow 



code of Johannesson et al. (2006) to calculate the lightcurves. It takes as input the physical 
parameters E- moi A*, An, R sw , p, e e , and e&, where An is the change in density at the wind 
termination shock and R sw is the radius of the termination shock. In addition to these 
parameters, we also fit for dust extinction Ay in the host frame. The multi-wavelength fit 
therefore directly provides the best fit physical parameters. Note that, in this model the 



definition of e e has been changed from the one used in Johannesson et al. (2006) to that of 



Panaitescu & Kumar (2001b) to allow for values of p < 2. 



3.1 Forward-reverse shocks 



Shao & Dai (2005) modelled the i?-band light curve alone using a forward-reverse shock 
model. In this model, for T <0.02 days the RS dominates the emission. The rise at T ~0.02 days 
coincides with the passage of the synchrotron frequency of the forward shock through R- 
band, after which the FS dominates the light curve. However, in the Swift V-b&nd, with 
a good sampling at T <0.02 days, the light curve decay is too shallow to be explained by 
standard reverse shock emission. Nevertheless we first attempted an forward-reverse shock 
model where the optical light curve at T <0.02 days is not entirely dominated by the RS, 
but also by the rising FS. This allowed us to have the steep decay of the RS, but still re- 
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produce the early optical decay. However, this placed constraints on the temporal profile of 
the modelled FS emission at the early period. To reproduce the observed optical and X-ray 
light curves, we had to have a wind driven density profile for the ambient medium, and a 



hard (p < 2) electron distribution (Bhattacharya 2001 Sagar et al. 2001 Misra et al. 2005 



Resmi & Bhattacharya 2008) undergoing fast cooling during this time period. We do not 



have to include dust extinction in the host frame to achieve this fit. Along with further 
constraints on the FS emission from the later part of the optical and X-ray light curves, 



and the Spitzer spectrum (Heng et al. 2008) at 2.3 days, we narrowed down the parameter 
space. In Figure [4| we present the best fit model (reduced chi squared, Xdof = 3-9 for DOF 
of 223) along with the optical, X-ray and radio data. We obtained a good fit in the spectral 
parameter space but the physical parameters inferred from these spectral fits turned out to 
be unrealistic (Table [3]) and based on that we had to rule out this model. 



3.2 Two-component outflow 

We attempt a two-component model, where different outflows dominate the early and late- 
time light curves. The most popular picture is of two co-aligned components with different 
opening ang les ( |Berger et al."|2003| [Racusin et al.|2008"j ). For GRB 050525 A, the first narrow 
jet dominates the emission at T < 0.02 days while the second (wide) jet gives rise to the 
.R-band re-brightening at T ~ 0.02 days as it decelerates around this epoch, and dominates 
thereafter. The late-time optical and X-ray light curves are seen to steepen gradually. We 
assume this to be due to the lateral expansion of the wide component and our best fit 
jet-break is at ~ 5 days. Due to the lack of dense sampling around this time, there are 
relatively larger uncertainties on this estimate. After T > 0.1 days the emission from the 
narrow component is required to decay faster so that the model does not over-predict the flux 
at later epochs. Even though various reasons could cause the quenching, the most natural 
assumption is that the narrow component is undergoing lateral expansion around that time. 
The narrow jet contribution is significant in X-rays for sufficiently long time and hence 
relatively good constraints on its tj can be obtained from the X-ray light curve alone. A 
constant density ambient medium profile produced the best fit. The model required a Milky 
Way type additional extinction in the host galaxy frame. 

In Figure [5j we present the best fit model along with the optical, X-ray and radio data. 
We obtain a somewhat high Xdof °f 4.5 (for DOF of 223). Such large values are due to the 
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small scale fluctuations in the light curves, especially the scatter around the re-brightening 
epoch in i?-band and the early flare-like feature in X-ray. Our best fit spectral and physical 
parameters are given in Table |4j Our fit parameters (especially the value of p, the host 
extinction and tj for the second jet) differ from the values reported in previous studies 



(Blustin et al. 2006 Delia Valle et al. 2006 Heng et al. 2008), because our numerical code 



uses a more sophisticated description than the multiple powerlaw fitting employed in these 
papers. Since the narrow component is relevant only for a short period of time, its spectral 
parameters are not well constrained. We are only able to obtain lower limits to the self 
absorption frequency (u a ) and cooling frequency (z/ c ). The inferred physical parameters are 



well within the range observed for other afterglows (Panaitescu & Kumar 2001a Resmi 



& Bhattacharya 2008 Cenko et al. 2010). The ambient medium density is towards the 



higher end of what is normally estimated in afterglow modelling, leading to a wider inferred 
collimation angle for the given value of tj. The estimated value of kinetic energy in both the 
jets are similar. The initial lorentz factor (r ) of the wide jet is ~ 10 and that of the narrow 
jet is > 60, suggesting the more relativistic narrow jet to be responsible for the prompt 
gamma-ray emission. The energy in the narrow jet is ~ 370 times that in the radiation 
observed by BAT (2.3 x 10 52 erg, ( Blustin et al.||2006 )). This could indicate that only around 
0.3% of total energy is converted to radiation during the prompt emission phase, which is 



compatible with the low efficiencies expected for internal shocks (Daigne & Mochkovitch 



1998). 



Since for both jets e e ^> e#, we estimated the expected inverse compton emission and 
found it to be significantly lower than the synchrotron component. 



The model agrees reasonably well with the VLA radio data in multiple frequencies, 
however the radio light curves are more or less flat for the entire duration of VLA observation. 
The 8.4 GHz and 4.8 GHz fluxes remain at the same level within error bars from 0.3 to 
30 days. This could be caused either by contamination from the associated supernova or 
by contribution from a nearby unresolved variable source (like an AGN). We note that a 
similar flattening has been observed in other radio afterglows ( Frail et al.|2004 ). Nevertheless, 
we have provided the best possible fits to the available radio data within our model. The 
steepening in the light curves corresponds to the epoch when the fireball becomes optically 
thin to synchrotron radiation. 

© RAS, MNRAS 000, [l]-?? 
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3.3 Wind driven bubble and termination shock 



The two-component model agrees well with the overall behaviour of the multi-wavelength 

light curves. However, for the early X-ray light curve (0.001 0.012 days), only the average 

behaviour is explained by this model. The X-ray light curve during this period displays mild 
flares. In an attempt to explain this finer feature, we opted for another model where the 
ambient density profile changes shape, such as in a wind bubble. The immediate neighbour- 
hood of the burst has a stratified density profile up to the termination shock R BW , where the 
afterglow shock encounters a density jump as well as a change in the density profile. When 
the forward shock enters the density enhancement it rapidly slows down and the co-moving 
density of the shocked matter is increased. There is therefore an increased number of radiat- 
ing electrons but their average energy is lower. This causes the afterglow temporal decay to 
be slower for v < v c but faster for v > v c because v c rapidly decreases when the shock enters 



the density jump (Johannesson 2006). The afterglow then settles into a powerlaw decay in 
the constant density medium. If v c remains between the optical and X-ray observing bands 
in this afterglow, this model can reproduce the optical bump observed in the light curve. 

Our best fit model is shown in Figure [6] where the wind termination shock is placed 
at ~ 0.09 pc (corresponding to the observer frame time T ~ 900 s), and has a density 
increase of a factor of ~ 12. The set of parameters used for the model are given in Table [5} 
The model parameters differ from those determined from the two-component model given in 
Table [4]because of the different modeling approach. The parameters are more like the narrow 
component in the two-component model, in agreement with the fact that it better fits the 
early afterglow data than the late data. The typical distance to a WR stars termination shock 
ranges from ~ 0.1 pc to a few hundred pc depending on the mass loss rate of the star and 
the interstellar medium surrounding it ( Eldridge et al.|[2006 van Marie et al.|20~0~6| pdridje 



2007b). The radius of the wind termination shock is much smaller than that expected for a 



typical Wolf-Rayet star, possibly indicating that the afterglow erupted from the rotational 



axis and the star was in a dense interstellar environment (Eldridge 2007a). 

While the model adequately describes the optical and infrared light curves and early 
X-ray light curve, it over-predicts the late-time X-ray light curve and under-predicts the 
late-time radio observations. A possible explanation for this difference is a change in the 
shock physics causing the cooling break frequency to decrease further, for example a small 
increase in the magnetic field strength of the shock. Further density fluctuations might also 
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Multi-wavelength modelling of GRB 050525 A afterglow 11 
explain the differences as numerical modelling have shown the wind of massive stars to be 

. We obtain a Xdof °f ~ 3.6 
(for a DOF of 223). For this model too, the scatter in i?-band around the rebrightening 
epoch contributes to the somewhat large value of \ 2 - 



more structured than we assume in our model ( Eldridge||2007a 



3.4 Supernova contribution 



Optical light curves of several low-redshift GRBs show a late-time red bump (Zeh et al. 



2004), often attributed to the emission from a supernova believed to be associated with the 



progenitor of the GRB (Heger et al. 2003). Around T ~ 5 days, a flattening of the i?-band 
light curve was observed for GRB 050525A (Delia Valle et al. |2006[ ) with a temporal profile 
similar to the 'SN bump' seen in other afterglows. The late-time (~ 36 days) VLT spectrum 



observed by Delia Valle et al. also showed similarities to the spectra of other GRB associated 



SNe. We compared the late-time i?-band light curve with the prototypical SN98bw (Galama 



et al. 1998), the first known example of a supernova associated with a GRB. We used the 



template optical light curves of SN98bw (Galama et al. 1998) k-corrected for the redshift 
of GRB 050525A. Best agreement with the observed i?-band data for T > 5 days were 
obtained after shifting the SN98bw template by —0.9 days and applying a temporal stretch 
factor of 0.8. We did not have to adjust the magnitude of the shifted SN98bw, implying the 
associated SN to be of similar brightness as SN98bw. The host galaxy i?-band magnitude 
was fixed at 25.1 ( Delia Valle et aL~||2006 ) . The SN component is shown as thin gray lines in 
Figures |3J [5] and |6j 



4 SUMMARY 

We have presented the optical afterglow of GRB 050525A in VRIJH photometric bands. 
Our data fill some gaps in the optical multi-wavelength light curves beyond 0.01 days and 
provide a better constraint on the optical decay index. Our IR observations, though confined 
to a narrow time bin, provide additional spectral constraints. The millimeter upper limits 
contributes to a better picture of the low frequency behaviour of the afterglow. 

We have undertaken a comprehensive multi-wavelength modelling of the afterglow, and 
tested various models against the data. The afterglow behaviour is too complex for a simple 
blast-wave model. We find that including emission from a possible reverse shock component 
is not sufficient to explain the afterglow evolution. Either the outflow is structured as a two- 
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Figure 1. The GRB 050525A field. The left image (at T ~700 s after the burst, V-band) was taken with the 0.2m telescope 
at the BOOTES1-B telescope located at the INTA-CEDEA station at El Arenosillo, Huelva, Spain. The right one, shown for 
comparison purposes, was taken at T ~2 days after the burst with the 1.5m telescope at IAA-CSIC Observatorio de Sierra 
Nevada (OSN) in the /-band. The green circle is 5'.'0 in radius and indicates the location of the afterglow in the two images. 

component jet or the ambient medium has a complex structure involving a variation in the 
density profile. Our two-component jet model is able to reproduce the overall behaviour of the 
afterglow, except the finer fluctuation in the early X-ray light curve. The wind termination 
shock model succeeds in explaining the early phase including the short time-scale features 
but deviates from the late-time data. The radio light curves are moderately well explained 
by both models. Afterglow modelling is necessary to unravel the nature of the outflow and 
the structure of the medium around the burst. The complexity of the observed light curves 
demands that the most realistic models are used. Dense sampling and long monitoring 
campaigns are also required in conjunction with afterglow modelling. 
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Figure 2. UBVRIJH, X-ray, radio and millimeter light curves of the afterglow of GRB 050525A. The dark and light coloured 
symbols correspond to the new data presented in this paper and the literature data respectively. The upper limits, including 
the millimeter ones presented in this paper, are represented as open triangles. The dash-dotted vertical line corresponds to 
T ~ 0.02 days (see text). The optical light curves are shifted arbitrarily for clarity. 
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Figure 5. Same as in figure [4] but for the best fit two-component outflow model. For the X-ray light curve, contribution from 
the narrow and wide components are shown as dashed and dash-dotted curves separately. 
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Figure 6. Same as in figure[5]but for the best fit wind termination shock model. 
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Table 2. Log of optical and NIR observations 
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Time since 
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: burst (Ty 
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Exp time 
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Table 3. Best fit (Xdof ~ ^*") s P ec tial and physical parameters of the forward-reverse shock model. No host extinction need 
to be used. However, the values of the derived physical parameters are unrealistic and hence we rule out this model. 

fa (Hz) v m (Hz) v c (Hz) f m (mjy) p q tj (day) 

X 

u 

I ~io 10 s.o±|| x io 15 7.osj:f;| x 10 11 15.82^ ^ i.eal ,;™ 1.2 ± 0.2 0.04l°;°* 

•d 
S 

fc -Eiso.52 (erg) A» e e e B 9j (°) E t ot,52 (erg) 



3t 2 9 71 x io-e i.6lg;°xl0-^( 3xl0 7 9g J 012 (2.5 ± 1.8) X 10? 0.11 



+ 0.03 
0.02 



2. L i _ x 10" 



(Hz) i/ m (Hz) (Hz) /,„ (mjy) p t x (day) 

2.1o 74 x 10 11 2.95l5'4 7 X 10 15 "'-t^'l 2.08 ±0.07 7.9lJ' g6 x 10" 



(l. ± o.2) x io -3 1073IH5 



Table 4. Best fit spectral and physical parameters of the two-component outflow model. The best fit A„(host) is 0.16 i 0.06 
mag (Milky Way type). The spectral parameters are quoted at the jet-break time of the respective components. Variabilities 
in the early optical and X-ray light curves result in a large value of Xdof ~ ^.5. 



u a (Hz) v m (Hz) v c (Hz) f m (mjy) p tj (day) t doc (day) 

I 
a 
o 
a 

a 

S (4.6+g^) x 10 10 (3.5lgg) X 10 10 (4.0 ± 3.0) x 10 16 2.2315'.? 2.6 ±0.08 5.0 ± 3.0 0.025 

o 

* £130,52 ( cr g) »0 (atom/cm 3 ) £ e e B Bj (°) Etot,52 (erg) T 

9-35lg° 3 354. S 2 ^l 6 0.126 ±0.02 (5.1 ± 3.4) x 10" B (24.81*;?) 0.43 ± 0.25 7.1 ±2.1 



v a (Hz) 



v m (Hz) 



(Hz) 



/m (mjy) 



tj (day) t doc (day) 



> 6.3 X 10 11 (I.oIqJ) X 10 11 > 2.6 x 10 



,+1.7 
-2.6 



2.37 ± 0.02 



. +0.02 
-0.01 



> 7.5 x 10" 



B 852.56 



Biao,52 (erg) n (atom/cc) 

1/4 



fij (°) Etot,53 (erg) 



(axiois) 



354.8 



°' 014 (^w) 1/4 6 - 1Xl0_6 (3^Trr) _5/4 < 3 - 26 "(Siw)' 



r 



> 60 
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Table 5. Parameters of the wind termination shock model. For explanation of para meters see l|J6hannesson et al, | | 2006 ) . 
Note that the definition of e e has been changed to allow for values of p < 2 as in (Panaitcscu & Kumar 2001b|. Values 
without uncertainties were fixed in the model and do not affect the results significantly. Dust extinction in the host frame is 
A v = 0.15 ± 0.05 (SMC type). The best fit model yields a Xdqf ~ 3 - 6 - 



Parameter Value 

Etct,52 (erg) 0.lt°;j> 8 

e } n < i.i 

A 5+ 10 

* U -°-0.45 

Ant > 12 

R\ w (pc) 0.09 ± 0.06 

p 1.85 ±0.15 

ee/10- 5 1±|. B 



t Fractional change in density at the termination shock 
+ Radius of the termination shock 
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